Shocks and slip systems: predictions from a theory of continuum dislocation dynamics 
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Using a recently developed continuum theory of dislocation dynamics, we derive three new predic- 
tions about plasticity and grain boundary formation in crystals. (1) There will be a residual stress 
jump across grain boundaries and plasticity-induced cell walls, which self-consistently acts to form 
the boundary by attracting neighboring dislocations; we derive the asymptotic late-time dynamics 
of the grain-boundary formation process. (2) At grain boundaries formed at high temperatures, 
there will be a cusp in the elastic energy density. (3) In early stages of plasticity, when only one 
type of dislocation is active (single-slip stage I plasticity) , cell walls will not form; instead we predict 
the formation of jump singularities in the dislocation density. 
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Dislocation motion governs the plastic deformation of 
crystals, as well as the formation and evolution of grain 
boundaries in polycrystals. Acharya and Roy recently 
proposed an evolution law for the density of dislocations 
in crystals, allowing both glide and climb (hence appro- 
priate at high temperatures). We subsequently rederived 
this law from the microscopic dynamics and a closure ap- 
proximation, proposed a modified law to suppress climb, 
and showed numerically that both laws developed shock 
singularities We argued that these shocks are the 
fundamental origin of the wall-like structures sponta- 
neously formed by dislocations — providing a continuum 
explanation of both the grain boundaries formed at high 
temperatures and the cell walls formed at low temper- 
atures in plastically deformed crystals. Wall formation 
thus emerges naturally from this continuum formulation. 
What other properties and predictions naturally emerge 
from this mesoscale perspective? Here we discuss three. 

(1) Residual stress. The continuum model predicts 
residual stresses associated with the formation of grain 
boundaries and cell walls. It makes detailed predictions 
about the asymptotic behavior of the stress and dislo- 
cation density near the walls, and about how the wall's 
associated stress allows it to draw in and absorb neigh- 
boring dislocations. 

(2) Cusps in energy density. The continuum model 
predicts that the elastic energy associated with the resid- 
ual stress is continuous across a grain boundary (despite 
the jumps in the individual stress components). The en- 
ergy density does have a cusp singularity. 

(3) Dislocation density jumps in single-slip. When ap- 
propriately restricted to forbid climb, the model predicts 
that materials with only one kind of dislocation will not 
form cell walls, explaining the lack of such structures in 
systems (discrete-dislocation simulations and experimen- 
tal stage I plasticity) where only one slip system is active. 
Instead, the model predicts that these systems will form 
a more subtle jump singularity in the dislocation density. 

We will derive these three predictions analytically by 



mapping the dislocation dynamics theory in one dimen- 
sion onto Burgers equation. We will present numeri- 
cal simulations to flesh out these predictions. Finally, 
we will consider how and whether these predictions are 
likely to depend upon the detailed structure of the contin- 
uum model, and how physical mechanisms ignored by the 
model will likely affect and modify these striking predic- 
tions. We begin by introducing our notation and Acharya 
and Roy's derivation of the equations of motion we study. 

A complete macroscopic description of the deformation 
u of a material is given by diUj = Pfj -\- (3fj, where Pfj 
represents the elastic, reversible distortion and the plastic 
distortion tensor f3fj describes the irreversible plastic de- 
formation. The plastic distortion is the result of the net 
density of dislocations ,_described by the Nye dislocation 
density tensor 
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which measures the net flux of dislocation a, tangent to 
t, with Burgers vector b, in the (coarse-grained) neigh- 
borhood of X. 

The time evolution of the plastic distortion tensor is 
given in terms of a function Jij = dtPfj, where eumJmj 
gives the net current for the conserved Burgers vector of 
the dislocation density tensor, dtPij = —eumdiJmj- The 
current from a single dislocation moving with velocity v is 
Jij — £iintibjVn6{^), and the net Peach-Koehler force on 



a dislocation driving its motion is /; 
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where cr is the local stress (due, for example, to the other 
dislocations). In our derivation of the equation of mo- 
tion 0, we allowed each dislocation to move indepen- 
dently, and then made a closure approximation to write 
J in terms of p. Roy and Acharya [l| got the same final 
result by simply assuming that all the dislocations move 
with the same velocity v, given by a constant (D{p)/2) 
times the force density J^i ~ —SimnPmc'^nc on the local 
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dislocations: 
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A physically natural choice for D{p) is proportional to 
an inverse density of dislocation lines (so that the force 
per dislocation drives the motion Q); we choose here to 
discuss the mathematically more convenient choice of a 
constant D. (Numerically, both give qualitatively similar 
evolution We controlled the microscopic mobility 

difference between glide and climb by introducing a 
mesoscale parameter A. By setting 

dtPjj ~ Jij = Jfj^ - -J 5i]Jkk^^ (3) 

at low temperatures A = 1 removed the trace of J enforc- 
ing volume conservation, and hence forbids climb, while 
at high temperatures A = allowed for equal mobilities 
for both glide and climb. 

We now turn to an analysis of the continuum equation 
in one dimension. Near a wall singularity (say, perpen- 
dicular to z) the dynamics is one-dimensional. The varia- 
tions of the stress, plastic strain, and dislocation densities 
parallel to the wall asymptotically become unimportant 
compared to the variations along z as one approaches the 
singularity. In one dimension the stress o"(x), generally 
given by a long-range integral over the neighboring dis- 
locations p(x'), can be written as a linear function of the 
local plastic distortion /3^: 

Cij = —CijkmPkm (4) 

where 

Cijkm = M (^ikSjm + SimSjk + SijSkni^ (5) 

(not quite equal to the elasticity tensor), and the Kro- 
necker delta is modified such that 6zz = 0. Hence the evo- 
lution law for simplifies dramatically. The (3^^ com- 
ponents do not evolve (except for A 7^ 0, which helps 
enforce volume conservation); the other components all 
evolve according to 

dtl3^ = -lid^S) dz (^(3^ - ^ (6) 

where we have rescaled the time to set D ^ 1. The elastic 
energy density £ in equation can be written in terms 
of the local plastic distortion tensor: 

^ = 2 ^^i'^^^l^ijf^km ~ '^iPxy + Pyx) 

+ M (C' + P^y") + Y^iPL + fyy? (7) 



We now specialize to the case of A = 0, where glide and 
climb are treated on an equal footing (applicable to grain 
boundary formation during polygonization at high tem- 
peratures, for example). In this case, equation El tells us 
that the individual components of (3^ are all independent 
of one another, slaves solely to the evolution of the total 
stress energy density E. By contracting equation |S1 with 
CijkmPkm then using expressiond the time evolution of 
the strain energy becomes 

dt£ + \{dz£f ^Q. (8) 

Equation|Slcan be cast into the famous Burgers equation 
by defining = H H [13 : 

dtT + TdzT ^ 0. (9) 

The scalar ^{z) is again the net Peach-Koehler force 
density on the local dislocation density p{z). Burgers 
equation is the archetype of hyperbolic partial differen- 
tial equations. Under Burgers equation J- will develop 
sharp jumps downward after a finite evolved time, corre- 
sponding to cusps in the energy density £ and leading to 
jumps in the components of (3^ (Fig[T|l. 
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FIG. 1: Cusps and jumps in one dimension. For the 

continuum dynamics equation allowing both glide and climb, 
the Peacli-Koehler force density T obeys Burgers equation, 
and hence develops sharp jumps (top). The individual com- 
ponents of the plastic distortion tensor (bottom), as well 
as the stress and strain tensors, evolve to also have sharp 
jumps at the walls; the dislocation density hence develops the 
5-function singularity associated with grain boundary forma- 
tion. These discontinuities in the residual stresses and strains, 
however, cancel out in the net elastic energy density, which is 
continuous with only a cusp at the walls (middle). 

At late times, the asymptotic late-time solutions to 
Burgers equation (between the singularities) are linear 
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functions of z whose slopes decay with time: 



t-to 



(10) 



The corresponding clastic energy density (continuous 
across the singularities) has asymptotic form given by 
integrating equation 1101 



^ 1 {z - zpf 
^ 2 {t- to) 



So 



(11) 



The individual components of the distortion tensor 
(apart from the three time-independent components f3zj) 
numerically take the form 



Z- Zq 



(12) 



which can be shown to be consistent with the evolution 
law for the energy density feauation llll) . so long as the 
coefficients aij and 7^ obey the relations 



C'ijkm^ij'^km — 1 j C'ijk^-fidij^kni — 0, 



(13) 



for i,j = x or j/; there are no restrictions on aiz and 7^^. 

There has been extensive work on discrete point dis- 
location simulations in two dimensions, where it appears 
necessary to include more than one slip system to form 
walls Ullllllllllllllllllllli- What does the 
continuum model predict for this case? Let us con- 
sider a 2D system (constant along z) of straight edge 
dislocations with Burgers vector along x, described by 
a single non-zero component pzx{x,y). Such a system 
has two non-zero components of the distortion tensor, 
Pzx = ~9x(3yx + dy(3^^. 

A simulation of this system with Gaussian random ini- 
tial /3^2. and (3^^, allowing both glide and climb, generates 
a series of walls of dislocations roughly parallel to the y- 
axis similar to those seen by Barts and Carlsson in 
their 2D study of single slip with both glide and climb. 
(See figure 13) 

To forbid climb in this case it is convenient |23| to 
simply choose (3^.^. = 0. Figure 12 shows the evolution of 
the distortion field for climb-free dynamics with a single 
slip system. In agreement with experiment and the dis- 
crete dislocation simulations |20l l21j . we observe no cell 
wall structures in single slip (which would correspond to 
jumps in in Fig [21). Instead, we find a network of 
surfaces exhibiting a striking new singularity: a cusp in 
the distortion tensor, corresponding to a jump in the dis- 
location density. We can understand this singularity an- 
alytically using our mapping to Burgers equation. With 
only one non-zero component of /3^, the energy density 
is proportional to the square of (iy^ (eqnEI). The energy 
density satisfies eqn|Sland forms cusps, so the distortion 
tensor must also form cusps. 




FIG. 2: Cusps formed with one slip system. Plastic dis- 
tortion tensor j3^^ formed by climb-free evolution of a Gaus- 
sian random initial state of edge dislocations pointing along z 
with Burgers vector along ±x. Notice that walls do not form 
with one slip system, only cusps in the distortion tensor; com- 
pare to [21! . 



Which of these new predictions of the mesoscale dis- 
location dynamics model can be trusted? Which seem 
(in retrospect) physically plausible, and which are likely 
artifacts of simplifications made in the modeling process, 
or just mathematical curiosities of this particular model? 

(1) Residual stress. The primary driving force for dis- 
location motion is stress. In cases where dislocations 
dynamically assemble into walls (polygonization at high 
temperatures, cell wall formation at low temperatures), 
it does seem natural in retrospect to expect that the 
walls will be associated with stress jumps designed to 
attract residual dislocations to the walls. Grain bound- 
aries formed from the melt when separate growing crys- 
tals touch should likely not be described by the contin- 
uum model |23l |. These residual stress jumps must be 
viewed as a fundamental prediction of the model. 

(2) Cusps in the energy density. If there are jumps in 
the stress at grain boundaries, surely it is natural that 
there be some singularity in the energy density. At late 
stages when the dislocations between grain boundaries 
have all been removed, a flat boundary can lower the sys- 
tem energy by moving into the region of higher energy 
density. If both glide and climb are allowed [l^l, and if 
dislocation mobility is unimpeded (by precipitates, im- 
purities, lattice pinning, or tangling) this traction will 
lead the boundaries to move until the energy density is 
continuous across the boundary. Hence, for mobile walls 
at high temperatures, it is natural to expect the energy 
density to be continuous, and have only cusp singulari- 
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FIG. 3: Continuum of walls. The dislocation density tensor 
pzx evolved allowing both glide and climb from a random 
initial state of edge dislocations along t ~ z with b ~ x. 
Notice that the dislocations arrange themselves into small- 
angle tilt boundaries at the lattice scale, but do not coarsen; 
compare to [l9||. 



ties. 

(3) Dislocation density jumps in single-slip. The for- 
mation of dislocation walls in the continuum theory can- 
not properly be called a prediction, since wall formation 
is well known experimentally. It is, however, a natural 
consequence of the hyperbolic form of the equations. The 
(hitherto unobserved) prediction of dislocation density 
jumps in early stages of plasticity when only one slip sys- 
tem is active is also a natural consequence of hyperbolic 
equations. Our continuum model is guaranteed to lower 
the net energy with time, reassuring us that the predicted 
dislocation density jumps are energetically favorable and 
satisfy all compatibility constraints. They could, how- 
ever, be smeared by pinning and inhomogeneities in real 
systems just as grain boundaries and cell walls are dis- 
torted by these effects. A smeared dislocation density 
jump may be more challenging to identify than a smeared 
wall of dislocations, perhaps explaining why these jumps 
have not yet been seen experimentally. 

We thus predict that residual stress is not due to inho- 
mogeneities or other infelicities in the formation process, 
but is an intrinsic component of the formation of grain 
boundaries and cell walls. We make concrete predictions 
about the nature and form of these stresses near grain 
boundaries, that should be testable in colloidal systems 
or using next-generation X-ray sources. We provide a 
mesoscale explanation for one of the key morphological 
distinctions between early (stage I) and later (stage III) 
plasticity. Finally, we predict a new dislocation density 



jump structure in plastically deformed systems with only 
one active slip system. 
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